Universal dissipation scaling for non-equilibrium turbulence 
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It is experimentally shown that the non-classical high Reynolds number energy dissipation be- 
haviour, Ce = eL/u^ = f(ReM)/ReL, observed during the decay of fractal square grid-generated 
turbulence is also manifested in decaying turbulence originating from various regular grids. For 
sufBciently high values of the global Reynolds numbers Rcm, /(R&m) ~ Rcm- 
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In recent papers describing the wind tunnel turbulence 
generated by fractal square grids [H 0] it was shown that 
the turbulent kinetic energy dissipation rate, e, at mod- 
erately high Reynolds numbers does not follow the ex- 
pected scaling eL/u^ = ~ const (where L is the longi- 
tudinal integral length-scale and u the streamwise r.m.s. 
velocity). Instead [ll Q found that = f{ReM)IReL 
during the turbulence decay where /{Rgm) is an increas- 
ing function of Rbm = UooM jv, a global Reynolds num- 
ber based on a length-scale M characteristic of the grid, 
and where Rcl = uLjv is a local, downstream position 
dependent, Reynolds number {v is the kinematic viscos- 
ity and Uoo is the inflow velocity). This behaviour is 
accompanied by a well-defined power-law energy spec- 
trum (with exponent close to Kolmogorov's -5/3) over a 
broad range of length-scales and is therefore caused by 
a physically different underlying phenomenon than the 
well-known low Reynolds number law ~ Rcj} ■ 

Evidence of such a non-classical behaviour is significant 
due to the central role the empirical law ~ const has 
on most, if not all, models and theories of both homo- 
geneous and inhomogeneous turbulence [^-Q. Clearly, 
one should expect the existing models to inadequately 
describe turbulent flows (or regions thereof) not obeying 
the Ce ~ const scaling and consequently fail in their pre- 
dictions of transport phenomena (energy transfer, dissi- 
pation, particle dispersion, scalar diffusion, etc.). Most 
importantly, it challenges our understanding of turbu- 
lence phenomena in general, nevertheless providing a 
starting point for its study as well. 

In this Letter we report results which show that this 
non-classical behaviour is in fact more general than pre- 
viously thought and is not exceptional to the very spe- 
cial class of inflow conditions defined by fractal square 
grids. Hence this non-classical behaviour is of general sci- 
entific and engineering significance and therefore of much 
greater importance. 

In the present experiments we compare turbulence 
generated by three different regular square-mesh grids 
(RG230, RG115 and RG60) with the turbulence gener- 
ated by the fractal square grid (FSG) of [ij (see Fig. 
[U and table H]) . Our aim is to investigate the origin for 
the non-classical dissipation behaviour of the FSGs. The 
dimensions of RG230 are purposefully similar to those 



of the largest square on the FSG. This allows a ceteris 
paribus comparison between RG230 and FSG in two re- 
spects: (i) comparable inflow Reynolds numbers Rcm 
for similar inflow velocities if M is taken to be the side- 
length of the largest square on the grid (see Fig. [IJ 
and (ii) comparable distance from the grid where the 
wakes of the RG230 bars meet and where the wakes of 
the FSG largest bars meet. Starting from any one of our 
grids, the turbulent kinetic energy increases as one moves 
downstream along the tunnel's centreline and reaches a 
peak at a streamwise distance oipeafc from the grid be- 
yond which the turbulence decays [H, d, 0|- This dis- 
tance Xpeak is closely related to the distance from the 
grid where the wakes (largest wakes in the case of FSG) 
meet. Indeed, ^] introduced the wake interaction length- 
scale = M^/to where to is the lateral thickness of the 
largest bars (see Fig. [T]) and showed that x^ak scales 
with a;* in the case of FSGs. Subsequently, |l[ showed 
that Xpeakjxi, took comparable values for RGs and SFGs, 
a point which the experiments reported in this Letter al- 
low us to confirm (see table |l|. The length-scales Xp^ak 
and X* turn out to be paramount for a meaningful com- 
parison between grids. 

There are of course important differences between the 
four grids used here, for example different values of block- 
age ratio a (ratio between the blocking area of the grid 
and the area of the tunnel's test section) and different 
values of x* (see table |l|. These differences cause differ- 
ences in various mean fiow and turbulence profiles across 
the tunnel section. However, they have no baring on our 
main finding that the outstanding behaviour previously 
found in FSG-generated turbulence is also present in tur- 
bulence generated by regular grids for a region whose ex- 
tent is determined by a;*. Beyond this region, in the one 
case (RG60) where we can reach sufficiently far beyond 
it as a result of the wind tunnel's test section being much 
longer than we find the classical behaviour Ce ~ const 
provided the Reynolds number is sufficiently high. 

The experimental apparatus described in [l| was re- 
peated for the present experiments with the length of the 
0.46m X 0.46m-wide test section shortened from « 4.5m 
to Ri 3.5m to match the extent of the longitudinal tra- 
verse mechanism. We also installed a grid at the entrance 
of the diffuser to maintain a slight overpressure across 
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TABLE I. Details of turbulence-generating grids; d is the lon- 
gitudinal thickness of the bars. 



Grid 




M to 


d 












(mm) (mm) 


(mm) 


(%) 


(m) 




RG230 


mono-planar 


230 20 


6 


17 


2.65 


0.63 


RG115 


mono-planar 


115 10 


3.2 


17 


1.32 


0.63 


RG60 


bi-planar 


60 10 


10 


32 


0.36 


:^ 0.4_^ 


FSG 


mono-planar 


237.7 19.2 


5 


25 


2.94 


0.45~ 



Taken from measurements of a very similar grid. 



the test section. All data are recorded with one- and 
two-component hot-wire anemometers operated at con- 
stant temperature. The main data are recorded with two 
in-house etched Pl-(10%)Rh single-wire {SW) sensors, 
SWl and 5*^1^2, having sensing lengths of Iw = 0.5mm 
and Iw = 0.2mm and wire diameters of dw = 2.5/im 
and dw = Ifim, respectively. A Dantec 55P51 cross-wire 
(XW) with Iw = 1.0mm and d^ = bfim is also used to 
record basic isotropy statistics. The spatial resolution 

of the measurements, quantified by l^/f] = {i^^/s) ^ 
is the Kolmogorov microscale; the isotropic estimate of 
dissipation e = 15i^{du/dx)'^ is used), is given in table [TT] 
for the furthermost up- and downstream locations and 
for the different inflow velocities. We repeated the elec- 
tronic tests to confirm that the maximum unattenuated 
frequency response of the SWs was at least kij = 1 {k is 
the wavenumber). The data acquisition and processing 
methodologies are also similar to those described in 
An exception is that we use, for simplicity, the classical 
Taylor's frozen field hypothesis to convert temporal into 
spatially varying signals, although we checked that this 
does not meaningfully affect the results. 

This Letter's new data are recorded along the centre- 
line in the lee of each of our four grids (Fig. [T]and tables 
HI and ITT)) . Data recorded between a grid and its corre- 
sponding Xpeak are excluded (see caption of table |lTl as 
we confine our study to decaying turbulence. In these 
decay regions, u/v (where v is the r.m.s. lateral veloc- 
ity) is typically between 1.2 and 1.1 and the ratio of the 
mean square of the lateral turbulence velocity derivative 
with respect to the streamwise coordinate x to the mean 
square of the streamwise turbulence velocity derivative 
with respect to x takes values between 1.4 and 1.6. Both 
ratios vary by less than 5% along the streamwise extent 
of our records. Note that Xpeak is about as long as half 
the wind tunnel's extent in the cases of RG230 and FSG 
(see table |T|. The RG60 was investigated in fT] where 
it was shown that for sufficiently high inlet velocities the 
dissipation followed a convincing « const during decay 
far downstream. We repeat those measurements using a 
higher resolution sensor {SW2) and include recordings 
much closer to the grid (table HI)). 

First, we compare the dissipation scalings of the de- 
caying turbulence originating from RG230 and FSG. The 



TABLE II. Overview of the experimental results. Xmin & 
Xmax are the first and last measurement locations correspond- 
ing to 0.48a;* & l.OQa;*, 0.64a;, & l.lQa;,, 0.61a;, & 2.38a;, and 
0.72a;, & 8.75a;, for FSG, RG230, RG115 and RG60, respec- 
tively. Probe SWl is used for the measurements of the first 
two grids and SW2 for the last two. 



Grid Symbol Uoo Reu u/Uoo{%) Rex 

(mS~^) (xlO^) X,mn I x„ 



Iw/rj 



FSG 



RG230 



♦ 

□ 

O 
★ 



RG115 ® 



RG60 



<1 

A 



15.0 
17.5 

5.0 
10.0 
15.0 
17.5 
20.0 

20.0 

10.0 
15.0 
20.0 



237 
277 

77 
153 
230 
268 
307 



9.7/5.0 



7.2/4.1 



385/249 4.8/3.0 
418/275 5.5/3.5 

180/140 1.8/1.3 
261/200 2.9/2.2 
326/258 3.9/3.0 
348/281 4.4/3.3 
385/300 4.9/3.7 



153 6.9/2.7 255/160 2.3/1.1 

40 177/96 2.8/0.6 

60 15/2.2 240/111 3.8/0.8 

80 290/135 4.7/1.0 




FIG. 1. Turbulence generating grids. 
FSGpJ], RG230, RG115 and RG60. 
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FIG. 2. Longitudinal energy density spectra Fn per 
wavenumber k of turbulence generated by RG230 for (black) 
Uoo = 20ms"^, a;/a;, = 0.64, (dark grey) Uoo = lOms"^, a;/a;, = 
0.64 and (light grey) Uoo = 5ms"^, a;/a;» = 1.19. 



3 



0.65 - ♦ ^ ■ 

\ 

♦♦\ 

0.6 - ^\ 

♦ 



^ 0.55 - 



2000 



4000 5000 







'J^ j X 



Re'^Jj^/Rex 



FIG. 3. Normalised energy dissipation Ce versus local 
Reynolds number Re^ of turbulence generated by FSG, 
RG230 & RG115 for different inflow Reynolds numbers Rem- 
The dashed lines follow o^ReJ^ for different Rcm- The Re\ 
values of the data in this plot range between 140 and 418. 



Reynolds numbers Re\ = uXjv (where A is the Taylor mi- 
croscale) at our measurement stations are given in table 
im and are all large enough for a significant separation 
to exist between the large, energy containing, eddies and 
the smallest dissipative eddies. Indeed, the scale sep- 
aration at the highest Reynolds number is L/ry « 460. 
The measured one-dimensional longitudinal energy spec- 
tra Fii exhibit clear power-laws over more than a decade 
with an exponent close to Kolmogorov's -5/3, at least for 
Reu > 2.3x10^ and Re\ > 250 (see Fig. [2] where we only 
plot RG230 spectra for brevity and clarity; FSG spec- 
tra can be found in [ij). However, both for RG230 and 
SFG, the cornerstone assumption of turbulence theory, 
Cg Ri const, does not hold in this region where the tur- 
bulence decays (between about 1.3m from the grid and 
the end of the test section) at these Reynolds numbers 
(see Fig. Instead, for any fixed Rcm, ~ Rej^ (as 
one moves along x) is a good qualitative approximation 
(in Fig. [3] each set of symbols corresponds to one Rcm 
and one grid, see table |lTl Rcl decreases as x increases). 
At the furthest downstream locations which correspond 
to the lowest Rcl values for each Rcm in Fig. [3l there 
is a slight departure from ~ R^l > probably due to far 
downstream test section confinement effects discussed in 
[l[. (In our records, L reaches a maximum value smaller 
than M/4 at Xmax for all grids.) Note that the well- 
known relation e = Ibvu^/X^ (e.g. [i^) and the definition 
of imply 15(L/A)2 = CeRcL and 15L/A = CeRex which 
means that Cg ~ Rej} is equivalent to ~ R^x^ and that 
such Ce behaviour implies L/A « const during decay. 

When, instead of keeping Rcm fixed and varying x, we 
keep X fixed and vary Rcm, we then find a very different 
dependence of on Reynolds number; asymptotically 
independent of it for both RG230 and FSG as Rem in- 



FIG. 4. Normalised energy dissipation Cs versus the Reynolds 
number ratio Re]^ /Re\ of turbulence generated by RG230 
and FSG for different inflow Reynolds numbers Rem- 



creases. If we keep with the usual expectation that 
is independent of v at high enough Rcm (which may be 
close to, but not exactly, true, see 0), then these two 
different dependencies on Reynolds number can be rec- 
onciled by 



Rcm 



Re J 



1/2 
M 



Rey 



(1) 



because u/Uoo and L/AI are independent oi Rcm to lead- 
ing order at high enough Reynolds numbers. Note that 
CgOcReMl Re L is equivalent to L/X 



1/2 

to Ce ocRcj^j /Rex- This equation is fairly well supported 
by our data both for FSG and RG230 at Ega/ > 2.3x10^^ 
(Fig. [?]) but with a grid-dependent constant of propor- 
tionality in ([T]). 

Equation ([T]) may appear to clash with the fact that Ce 
is approximately independent of both x and Re m in the 
case of RG60 at distances greater than about 1.5m from 
that grid in a wind tunnel test section of exact same 
width as the present one (see Fig. 7 in 111]). This is 
a distance greater than about Ax^ from the grid because 
x^ « 0.36ni for RG60. However, ([T]) has so far been estab- 
lished for decaying turbulence originating from RG230 
and FSG up to downstream distances of less than about 
l.5xt {xt takes much greater values for these grids, see 
table HI . It is therefore reasonable to investigate whether 

1 /2 

^ and its equivalent relation L/X ~ R^m hold at dis- 
tances below a few multiples of from the RG60 grid. 
In Fig. [5]we plot L/X as a function of the local Reynolds 
number Rex for RG60 at different levels of Rcm- We 
find that L/A « const in the region between 0.72a;* and 
2x* (where Rex takes the largest values) and that L/X 
and Rex decay in exact proportion to each other (i.e. 
L/X ~ Rex which is equivalent to Ce = const) at further 



1/2 

Rcj^j and therefore 
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FIG. 5. I//A versus the local Reynolds number Rex of tur- 
bulence generated by the RG60 for different Rem and by 
RG115 & RG230 for the same Rcm- The dashed line follows 
Ce/15Rex with Ce: = 0.92. 
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FIG. 6. Normalised energy dissipation versus the Reynolds 
number ratio Rbm/Rsl of RG60-generated turbulence. 



downstream distances, i.e. where x > 2a;*. The region be- 
tween 0.72x» and corresponds to the 10 highest Re\ 
data points in Fig. [S]for each Rcm- The x- independent 
(therefore i?eA-independent) value of L/A in this region 
is an increasing function of Rcm as implied by dD). Such 
L/A behaviour was previously reported only for FSGs 
P, Q and is now shown to be more general. Replotting 
the RG60 data so as to directly compare with ([T]), we ob- 
tain Fig. [HI Equation ([T]) is a fairly good representation 
of the data up to Rbm/Rgl = 50, i.e. in the turbulent 
decay region closest to the grid up to a; « 2x». At stream- 
wise distances larger than 2a;,^ where Rcm/Rsl is larger 
than 50, becomes approximately independent of both 
X and ReM as already observed in earlier studies (e.g. 

Our measurements of decaying turbulence originating 



from RG115 were designed for a direct comparison with 
RG230 at equal a and Rcm = 1.53x10^ but different mesh 
size AI. The data obtained from these measurements are 
reported in Figs. [3] and [S] and show that L/A and 
take effectively same values for the two grids and that 
these values are consistent with Ce = f{ReM)lR&L and 
constant L/X = ^/ f{ReM)/^5 in the ranges of x probed. 
However, Rbm is too low for ([J) to hold. 

The present data and those of [J, [2| conspire to form 
the conclusion that, irrespective of the turbulence gener- 
ating grid (Fig. [T|) and for high enough Rbm, 



(2) 



and equivalently L/A Ri \J Cii?eA//15 are acceptable ap- 
proximations in the non-equilibrium decay region Xpeak < 
X < Xe where for RG60 and Ci is a dimension- 

less constant which only depends on inlet/boundary ge- 
ometry (type of fractal/regular grid, a, etc). We might 
expect Xe to scale with a;* for other grids as well, and the 
equilibrium dissipation scaling e = C2u'^/L (where C2 
is an inlet/boundary geometry-dependent dimensionless 
constant, see [1, to be recovered at x > Xg for other 
grids too. However, our RG115, RG230 and FSG data 
and those of [l|, Q do not allow us to test these expecta- 
tions, nor do they allow us to explore how Xe/x* may de- 
pend on inlet/boundary conditions. RG230 and FSG, in 
particular, act as magnifying lenses which make the non- 
equilibrium region to be longer than the entire tunnel 
section's length. Equations ([T]) and ([2]), and more gener- 
ally Ce = f{ReM)/ReL which also covers lower values of 
Rcm, are approximately true in the non-equilibrium re- 
gion irrespective of flow/turbulence profile details which 
differ from grid to grid. The FSGs are magnifying lenses 
with added capabilities for tailoring flow and turbulence 
profiles which go beyond variations in a. 

Finally, it is important to stress that the energy 
spectrum has a well-defined power-law shape over nearly 
two decades with exponent close to -5/3 at the closest 
point to the grid that we sampled in the non-equilibrium 
region (Fig. [2]). This power- law region becomes progres- 
sively narrower with an exponent progressively further 
away from -5/3 as x increases. In the equilibrium region 
of RG60 where e ~ u'^/L, the energy spectrum is far 
from Kolmogorov-shaped. This may just be a conse- 
quence of the low Reynolds numbers in the equilibrium 
region of our RG60 runs. But it is remarkable that a 
near-Kolmogorov power-law shaped energy spectrum 
does in fact appear well before the turbulence has had 
the time to reach equilibrium. A similar observation was 
made in [lo| where near-Kolmogorov power-law energy 
spectra were reported in a cylinder wake within one 
cylinder diameter from the cylinder. 
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